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We report on progress and future plans for calculating kaon weak matrix elements for e'/e using staggered 
fermions. 



The first goal of this project is to provide a 
check of the results (e'/e < 0) obtained using 
quenched Domain Wall fermions (DWF) by the 
CP-PACS § and RBC § collaborations. The 
second goal is to extend the work to dynamical 
simulations, since quenching is likely one of the 
main systematic errors in present calculations j| . 
The third goal is either to find a window for new 
physics or to confirm the standard model, when 
our numerical results are compared with the ob- 
served values. 

Staggered fermions are appropriate for this 
purpose, because they preserve enough chiral 
symmetry to prevent composite operators from 
mixing with wrong chirality operators and to pro- 
tect quark mass from additive renormalization, 
which is essential to calculate e'/e. They have 
the advantages over DWF of requiring less CPU 
time and dynamical simulations are already pos- 
sible with relatively light quark masses. By con- 
struction, staggered fermions carry four degener- 
ate flavors (also called "tastes"). They allow large 
flavor changing quark-gluon interactions, which 
makes it impractical to non-perturbatively deter- 
mine matching constants for operators of our in- 
terest. Another disadvantage is that the unim- 
proved action and operators receive large pertur- 
bative corrections at one loop and have large scal- 
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ing violations of order a 2 . Both of these disadvan- 
tages can be alleviated by improving staggered 
fermions using smeared links. Flavor symmetry 
breaking in the pion spectrum is significantly re- 
duced with such links |^,||. 

Our project has progressed through the follow- 
ing steps. 

1. We calculated the current-current diagrams 
for the gauge-invariant unimproved staggered 
fermions at one loop ||, which provided a com- 
plete set of perturbative matching formula for 
e'/e, combined with existing results for penguin 
diagrams 0. 

2. We performed a numerical study on e'/e us- 
ing the Columbia QCDSP supercomputer ||. Us- 
ing the matching formula given in HQ], we con- 
structed fully one-loop matched gauge-invariant 
operators. As expected, we found large perturba- 
tive corrections for Bq, so that we cannot quote 
quantitative results for this although the statisti- 
cal uncertainty is under control. An unexpected 
result was that different quenched transcriptions 
of the continuum operators on the lattice (pro- 
posed by Golterman and Pallante M) lead to sub- 
stantially different values for Bq ||. This indi- 
cates a large quenching uncertainty and deserves 
further study. 

3. The goal was to find an improvement scheme 
which can reduce perturbative correction down to 
10% or smaller. To achieve this goal, we calcu- 
lated, explicitly, one loop matching factors for a 
variety of improved staggered actions and opera- 



tors §: 1) Fat 7, 2) Fat7+Lepage, 3) HYP, and 
4) Asqtad-like (Fat7+Lepage+Naik) actions. We 
observed that all the above improvement schemes 
significantly reduce the size of the matching co- 
efficients. After a higher level of mean-field im- 
provement, the HYP and Fat7 links lead to the 
smallest one- loop corrections |J. Since the HYP 
action reduces the non-perturbative flavor sym- 
metry breaking more efficiently in the pion spec- 
trum, we adopted the HYP scheme in our numer- 
ical study. 

4. We studied further on the HYP link (lOj. 
The HYP links possess some universal properties, 
which are summarized in the following 5 theo- 
rems. 

Theorem 1 (SU(3) Projection) 

Any fat link can be expanded in powers of gauge 
fields (A^). 

B, = B^+Bf+Bf)+... 
BW =0{A n ) 

1. The linear term, B^f 1 is invariant under 
SU(3) projection. 

(2) 

2. The quadratic term, is antisymmetric 
in gauge fields. 

Theorem 2 (Triviality of renormalization) 

1. At one loop level, only the B^ term 
contributes to the renormalization of the 
gauge-invariant staggered fermion opera- 
tors. 

2. At one loop level, the contribution from 
4 n) for any n > 2 vanishes. 

3. At one loop level, the renormalization of the 
gauge-invariant staggered operators can be 
done by simply replacing the propagator of 
the A u field by that of the B^ held. 

This theorem is true, regardless of details of the 
smearing transformation. 

Theorem 3 (Multiple SU(3) projections) 

1. The linear gauge field term Bjp in the per- 
turbative expansion is universal. 



2. In general, the quadratic terms may be dif- 
ferent from one another. But all of them 
are antisymmetric in gauge Relds. 

3. This theorem is true, regardless of the de- 
tails of smearing. 

Theorem 4 (Uniqueness) 

If we impose the perturbative improvement condi- 
tion of removing the flavor changing interactions 
on the HYP action, the HYP link satishes the 
following: 

1. The linear term B^ in perturbative expan- 
sion is identical to that of the SU(3) pro- 
jected Fat 7 links. 

(2) 

2. The quadratic term B n is antisymmetric 
in gauge helds. 



Theorem 5 (Equivalence at one loop) 

If we impose the perturbative improvement con- 
dition to remove the flavor changing interactions, 
at one loop level, 

1. the renormalization of the gauge invariant 
staggered operators is identical between the 
HYP staggered action and those improved 
staggered actions made of the SU(3) pro- 
jected Fat 7 links, 

2. the contribution to the one-loop renormal- 
ization can be obtained by simply replacing 
the propagator of A^ by that of B u . 

The first two theorems were used in |ll|] and Q , 
although they did not present their derivations. 
For derivations of all five theorems and further 
details, see [|l0|. As a result of these theorems, 
we can prove that for each Feynman diagram, 

II C/«t II < II C tMn || . (1) 

Here, Cf a t {CtMn) represents perturbative cor- 
rections to gauge-invariant staggered operators 
constructed using SU(3) projected fat links (thin 
links). This inequality is not valid for those fat 
links without SU(3) projection. Hence, this lead 



Operators [P x P][P x P] n 



NAIVE 2 x (111.3 - 2C N ) 
HYP 2 x (8.25 - 2C H ) 



Operators [S x P][S x P] n 



NAIVE 2 x (95.6 - 6CW) 
HYP 2 x (14.6 - 6C H ) 



Table 1 

One-loop correction to (O3)//. 

to a conclusion that we may view the SU(3) pro- 
jection of fat links as a tool of tadpole improve- 
ment for the staggered fermion doublers |L0| . We 
also present alternative choices of constructing fat 
links to improve staggered fermions in The 
above five theorems make the perturbative calcu- 
lation simpler for the HYP scheme, because one 
can perform the calculation merely by replacing 
the thin link propagator with that of the HYP 
links. This simplicity is extensively used in calcu- 
lating the renormalization constants of the four- 
fermion operators in the next stage. 

5. We calculate the current-current diagrams to 
obtain perturbative matching coefficients for the 
staggered four-fermion operators constructed us- 
ing the HYP links. Especially, we are interested 
in the (Oz)n operator (we use the same notation 
as in g), because this receives large perturba- 
tive corrections (« 1) in the case of unimproved 
staggered fermions. 

(0 3 ) u = 2([P xP][PxP]-[Sx P][S x P]) n 

In Tables |l| and |^, we present values of one-loop 
corrections to (03)// constructed using both the 
HYP link and unimproved thin link (denoted as 
NAIVE). Here, C N = 13.159 and C H = 0.7709, 
which correspond to tadpole improvement contri- 
butions. The results shows that, by choosing the 
HYP scheme, the perturbative corrections are re- 
duced to w 10% level. In the case of the HYP 
scheme, note that the size of one-loop correction 
is already under control even without tadpole im- 
provement. For further details, refer to fl3|| . 

6. At present, we are performing a numeri- 
cal study using the Columbia QCDSP supercom- 
puter. We calculate weak matrix elements for e'/e 
using the HYP links at quenched j3 = 6.0 on the 
16 3 x 64 lattice. So far, we have completed mea- 
surements on about 80 gauge configurations. 

We plan to complete the one-loop calculations 
with HYP fermions, including the penguin dia- 



Table 2 

One-loop correction to (O^u. 

grams [ju|. This will allow us to match the com- 
plete set of AI =1/2 and 3/2 four-fermion oper- 
ators constructed using HYP links to the corre- 
sponding continuum operators in, say, the NDR 
scheme. 

We have also generalized the definition of the 
optimal matching scale, q* H], to the case of 
matching factors for operators with non- vanishing 
anomalous dimensions jL3|. Ultimately, we wish 
to extend the calculation to partially quenched 
and unquenched QCD. 

We thank N. Christ, C. Jung, C. Kim, G. Liu, 
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